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ABSTRACT 

In the local universe, 10% of massive elliptical galaxies are observed to exhibit a peculiar property: a substantial 
excess of ultraviolet emission than what is expected from their old, red stellar populations. Several origins for 
this ultraviolet excess (UVX) have been proposed including a population of hot young stars and a population of 
old, blue horizontal branch or extended horizontal branch (BHB or EHB) stars that have undergone substantial 
mass loss from their outer atmospheres. We explore the radial distribution of UVX in a selection of 49 nearby 
E/SO-type galaxies by measuring their extended photometry in the UV through mid-infrared (mid-IR) with the 
Galaxy Evolution Explorer {GALEX), the Sloan Digital Sky Survey, and the Wide-field Infrared Survey Explorer 
(WISE). We compare UV/optical and UV/mid-IR colors with the Flexible Stellar Population Synthesis models, 
which allow for the inclusion of EHB stars. We find that combined WISE mid-IR and GALEX UV colors are 
more effective in distinguishing models than optical colors, and that the UV/mid-IR combination is sensitive to 
the EHB fraction. There are strong color gradients, with the outer radii bluer than the inner half-light radii by 
~1 mag. This color difference is easily accounted for with an increase in the BHB fraction of 0.25 with radius. 
We estimated that the average ages for the inner and outer radii are 7.0 ± 0.3 Gyr, and 6.2 ± 0.2 Gyr, respectively, 
with the implication that the outer regions are likely to have formed ~ 1 Gyr after the inner regions. Additionally, 
we find that metallicity gradients are likely not a significant factor in the color difference. The separation of color 
between the inner and outer regions, which agrees with a specific stellar population difference (e.g., higher EHB 
populations), and the ~0.5-2 Gyr age difference suggests multi-stage formation. Our results are best explained by 
inside-out formation: rapid star formation within the core at early epochs (>4 Gyr ago) and at least one later stage 
starburst event coinciding with z ~ 1 . 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation - infrared: galaxies - 
ultraviolet: galaxies 


1. INTRODUCTION 

Early-type galaxies (ETGs; ellipticals and lenticulars) are the 
oldest class of galaxy, and the final stage of galaxies formed 
at early epochs. Through multiple epochs they evolved into the 
most massive nearby galaxies, comprising at least 50% of the 
local stellar mass (Bell et al. 2003). The most massive of these 
are older and more metal rich (Nelan et al. 2005). Additionally, 
modified hierarchical simulations show that the most massive 
ellipticals have relatively shorter formation timescales, more 
progenitors, and later assembly times (De Lucia et al. 2006). 
Down-sizing scenarios (Cowie et al. 1996) predict that the stars 
formed at earlier epochs over a shorter time period. 

Since the majority of the luminosity for ETGs arises from 
the red part of the spectrum, they have been characterized as 
“red and dead.” Still, approximately 80% of nearby ETGs are 
detected in the ultraviolet (UV, i.e., GALEX/NUV; Schawinski 
et al. 2007). The primary source of the UV emission is cur- 
rently unclear. Stellar synthesis models have yet to converge on 
the combination of stellar properties that lead to a significant 
number of post-main-sequence stars that are long-term produc- 
ers of UV, creating a UV excess (UVX). In ETGs, the UVX 
appears as a sudden increase in the flux in spectral energy dis- 
tributions blueward of 2500 A, and was originally named the 


“UV upturn” (e.g., Code & Welch 1979; O’Connell et al. 1992; 
Brown et al. 1995). The origin of the UVX has been debated 
for several decades, with possible sources being recent star for- 
mation, metal-poor horizontal branch (HB) stars >10 Gyr old, 
younger metal-rich HB stars, or a combination of these (see 
O’Connell 1999, for a review). 

Methods for discerning the properties between red and 
dead ETGs versus galaxies with more recent star formation 
include: (1) photometric narrow- and wide -band properties, 
e.g., Ha, color-magnitude relations; (2) structural properties, 
such as detection of disk remnants or merger signatures; and 
(3) spectroscopic detection of age-dependent lines, such as Mg, 
[Oil], and H/L The presence of recent star formation deduced 
through ionized gas emission in ETGs has been examined for 
decades. There is little question that some ETGs have signatures 
of recent star formation in the their nuclei and even disk 
remnants, but detecting this emission requires sensitive, high- 
resolution imaging or spectroscopy, especially at the outer radii. 

The general consensus in the literature is that the source of 
UVX is likely older, extreme stars (>2 Gyr) because they create 
the temperatures needed based on UV slope. However, evidence 
for recent star formation in ETGs has also been observed (e.g., 
Salim et al. 2005; Yi et al. 2005; Schawinski et al. 2007; Kaviraj 
et al. 2007, 2009). The presence of trace star formation in ETGs 
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is generally correlated with environment, is mostly detected in 
low-density environments (e.g., Sarzi et al. 2007; Shapiro et al. 
2010 , and references therein), and is organized into disk-like 
and ring configurations. 

There is no clear spectral distinction between the post-main- 
sequence or low-mass stellar populations that can produce a 
large amount of UV emission. Blue post-main-sequence stel- 
lar populations, such as blue horizontal branch stars (BHBs) 
with extreme temperatures and low-mass stellar envelopes, are 
promising candidates. At higher temperatures (T e g > 14,000 K) 
BHBs are called extended (sometimes “extreme”) HBs (EHBs; 
Code & Welch 1979; Burstein et al. 1988; Greggio & Renzini 
1990; Ferguson & Davidsen 1993; Brown et al. 1995, 2003; 
Atlee et al. 2009; Conroy & Gunn 2010; Yi et al. 2011). 
Dorman et al. (1993) found that the driving factor for increasing 
the temperatures of BHBs is the mass of the envelope in the zero- 
age HB phase, where a smaller envelope allows for higher tem- 
peratures; higher metallicities will result in a bimodal distribu- 
tion of temperatures for envelope masses of 0.05-0.15 M G . The 
morphology of the HB in the color-magnitude diagram is crit- 
ically determined by the Reimers mass-loss parameter ( 17 ; e.g., 
Percival & Salaris 2011), which describes the efficiency of the 
mass loss of the outer envelope in the HB phase. The complex- 
ities of the BHB phase are not easily isolated with broad-band 
photometry or low-resolution spectroscopy. It is difficult to dis- 
tinguish the BHB stars from the upper main-sequence stars using 
line strengths. These complications are in addition to the famil- 
iar age-metallicity degeneracy (Lee et al. 2000; Trager et al. 
2005; Conroy & Gunn 2010). 

A metal-rich star would produce UVX through significant 
mass loss in the red giant branch phase: a resulting lower mass 
HB (asymptotic giant branch, AGB manque) star could sustain 
20,000 K for ~10 Myr (e.g., O’Connell 1999). Low- 
metallicity stars may also play a part since the oldest populations 
are thought to be metal poor; however, models suggest that 
the UVX is not present until >15 Gyr for Z < 0.5 Z Q (Yi 
et al. 1998). Observations show a lack of a UV upturn in metal- 
poor galactic clusters as compared with giant ellipticals (Yi 
et al. 1998 and references therein). Additionally, models that 
include UV upturn show that solar, or higher metallicities allow 
UVX at younger ages (<5 Gyr) compared with metal-poor 
populations (> 10 Gyr; see Figure 9 in Yi et al. 1998). However, 
this is strongly dependent on the mass-loss efficiency. It would 
require a significant amount (>20% of HB stars) to produce the 
observed FUV in galaxies (O’Connell 1999), making the ages 
and lifetime of EHB phases an important factor. In globular 
clusters, up to 30% of the HB stars appear to be in an EHB 
phase (Kaluzny & Udalski 1992; Liebert et al. 1994). 

Binary stars provide a promising alternative since their 
environments can host the conditions for rapid loss of the 
hydrogen envelope, forming hot helium-burning stars. There 
are no metallicity requirements for producing UVX in binaries. 
Observed candidates for such binaries include hot subdwarf 
(sdBs) stars in globular clusters (see review by Heber 2009). 
Han et al. (2007) model the UVX produced by binary sdBs in 
which stellar populations as young as 1 Gyr can reach these 
extreme temperatures and reproduce the colors seen for ETGs. 

Knowledge about how massive ellipticals evolve to their 
present state is important for understanding galaxy evolution. 
The theories need to integrate many factors including the mor- 
phologies and age-metallicity differences observed between 
more and less massive ETGs. Mergers play an important role in 
explaining these factors. For example, dissipational (gas-rich) 


mergers would seed new star formation, while a dissipationless 
(gas-poor) merger would add mass via combining older stellar 
populations. An older stellar population in the outer halo may 
also be explained by outside-in cessation, where the star for- 
mation ceases first in the disk (e.g., Pipino & Matteucci 2004; 
Pipino et al. 2006). In an inside-out cessation, a core forms 
through very rapid starbursts at high redshifts, and through ma- 
jor multi-stage growth at later epochs, the outer regions form 
via minor mergers and accretion of hot gas from the immediate 
environment (e.g., Daddi et al. 2005; Nelson et al. 2012; Saracco 
et al. 2012 ). 

The cores of nearby ETGs (within a half-light or effective 
radius) have been studied in much greater detail than the 
extended regions. The most recent and comprehensive study, 
the SAURON survey (de Zeeuw et al. 2002), uses the SAURON 
integral-field spectrograph to survey 72 nearby galaxies; 48 are 
ETGs. Among the sample, a kinematic separation is identified 
that separates the ETGs into subpopulations of slow and fast 
rotators based on their projected angular momentum (Emsellem 
et al. 2007). The slow rotators tend to be less compact and 
composed of older stellar populations (> 8 Gyr; McDermid et al. 
2006). Additionally, a slight trend with Mg 11 and FUV — NUV 
color is found for the slow rotators. The fast rotators have ages 
of <5 Gyr. Shapiro et al. (2010) use mid-infrared (mid-IR) 
colors with Spitzer to trace low-level star formation through 
polycyclic aromatic hydrocarbon (PAH) emission in the sample. 
Most significantly, they find that the slow rotators have no sign 
of PAH emission and are quiescent over an ~8 Gyr period. 

In summary, much literature has been written to address the 
origin of UV emission in ETGs, and how they may have evolved 
into their present form. Resolved spectral analysis on ETGs has 
mostly focused on scales slightly larger than the effective radius, 
because of practicality. For example, Jeong et al. (2012) use 
resolved spectroscopy to test the Burstein relation on galaxies 
with extended UV. However, the Mg lines used are within the 
effective radii. Another example is the work of Carter et al. 
(2011), who compare the UV and near-IR colors for the total 
flux and central regions of ETGs to determine if the FUV excess 
is related to velocity dispersion, metallicity, or abundance ratio. 
However, their radial profiles do not extend beyond ~30" from 
the galactocentric radius of each galaxy. These studies and 
similar ones do not use stellar population modeling to study the 
central and outer regions simultaneously to the spatial extent 
and resolution that is possible with photometry. In this paper, 
we conduct such an analysis as the first step in an expanded, and 
detailed, series of papers. 

We use mid-IR imaging from the Wide-field Infrared Survey 
Explorer {WISE) to explore the UV to mid-IR color-space 
of ETGs with UVX. The two main objectives in our study 
are to (1) determine a how photometry between the UV and 
mid-IR can be used to interpret the likely stellar populations 
contributing to the UV emission and (2) use radial information 
from these results to interpret how these galaxies may have 
evolved. We do this over the total flux of the galaxies in 
our sample to analyze the color differences between the inner 
and outer regions. WISE provides the all-sky capabilities for 
a complete study. In Section 2, we discuss our selection and 
the detailed extended photometry. In Section 3 we describe the 
stellar population synthesis models and parameters used for our 
analysis. We explore the color gradients, spatial distribution, 
and possible source of the UVX in our sample in Section 4. 
We conclude with a discussion on the physical origin and the 
possible formation scenarios that may explain our results in 
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Section 5. Where necessary, we adopt a ACDM cosmology, 
(An, A\, Ho) = (0.27, 0.73, 71 kins” 1 MpcT 1 ). 

2. SAMPLE SELECTION 

We selected ETGs from the GALEX - Ultraviolet Atlas of 
Nearby Galaxies (Gil de Paz et al. 2007). This atlas includes 
galaxies from the GALEX Nearby Galaxy Survey and galaxies 
within that field and other fields with similar or greater depth. 
They select galaxies based on optical parameters from de 
Vaucouleurs et al. (1991) with /x B = 25 mag arcsec -2 isophote 
diameters larger than 1'; they also use the de Vaucouleurs 
et al. (1991) galaxy types to label the morphologies. The atlas 
contains a total of 1034 galaxies where 893 have both FUV and 
NUV detections. For our sample, we selected by morphological 
type E or SO, resulting in a sample of 125 galaxies all at z < 0.06. 
We then applied three further criteria. First, since we selected 
the sources during the execution of the WISE survey, not all 
sources had full-depth WISE data at the time the analysis was 
performed. We therefore removed the sources without full-depth 
WISE data. This amounted to a random cull of the parent sample. 
Second, we removed sources that did not have coverage in 
both Galaxy Evolution Explorer {GALEX) channels. Third, we 
removed sources that lay near bright foreground objects, or had 
obvious artifacts in the data. This resulted in a final sample of 49 
reliable objects at a mean redshift of 0.02. Since the selections 
we applied to the Gil de Paz et al. sample are effectively random, 
our final sample should be representative of the local UV- and 
mid-IR-bright ETG population as a whole. 

2.1. Data 

We cross-matched our sample with three catalogs: the Sloan 
Digital Sky Survey (DR7 SDSS; York et al. 2000, u, g, r, i, and 
Z filters) the Two Micron All Sky Survey (2MASS, A -hand; 
Skrutskie et al. 2006), and WISE (3.4, 4.6, 12, and 22 /xm 
filters; Wright et al. 2010). This provides photometry in up to 
14 bands, over 0.15-22 /xm. The WISE point-spread functions 
for co-added images in 3.4, 4.6, 12, and 22 /xm are 8.3, 9.1, 9.5, 
and 16.8 arcsec, respectively. The sample comprises 25 SO and 
24 E galaxies with a mean 3.4 /xm radius of R tot = 55 kpc. 

To obtain integrated magnitudes, we proceed as follows. 
For 2MASS, we obtain data from the extended source catalog 
(2MASS XSC; Jarrett et al. 2003). For WISE we obtain 
magnitudes from the extended source catalog project {WISE 
XS), which is described in Jarrett et al. (2013) and briefly here. 
For each galaxy, the total flux is integrated over an elliptical 
surface determined by the convergence of the surface brightness 
profile and the background. The background was estimated by 
an outer annulus, the same minor/major-axis ratio and position 
angle {b/a and PA.) as the elliptical aperture. Each image was 
visually inspected to remove foreground stars. 

For GALEX and SDSS, we follow the methods in Neill et al. 
(2011). In the same method as Neill et al. (2011), we obtain 
archival GALEX images to coadd for the deepest possible image. 
We used both GALEX imaging filters (Martin et al. 2005): FUV 
(Left = 1539 A) and NUV (L eff = 2316A). For the SDSS 
data, Neill et al. (2009, 2011) found the SDSS catalog data to be 
inaccurate due to the galaxies extending over multiple strips and 
consequently having been broken into subregions. The authors 
developed a coadd and mosaic technique to obtain the proper 
integrated fluxes over multiple frames. We use these fluxes, and 
further correct for Galactic extinction using the dust maps of 
Schlegel et al. (1998). 


Table 1 

Early-type Galaxy Sample Properties 


Object 

R.A. 

Decl. 

fo/a a 

«tot b 

(kpc) 

«h b 

(kpc) 

cz c 

(km s" 1 ) 

Type d 

NGC 155 

00:34:40.0 

-10:45:58.7 

0.87 

75 

8.1 

6210 

SO 

NGC 163 

00:35:59.7 

-10:07:17.5 

0.93 

60 

7.7 

5982 

E 

IC 1698 

01:25:22.0 

14:50:19.7 

0.41 

27 

8.7 

6503 

SO 

IC 1700 

01:25:24.5 

14:51:52.8 

0.81 

58 

8.6 

6356 

E 

UGC 1040 

01:27:35.8 

-01:06:18.3 

0.25 

32 

6.3 

4489 

SO-a 

NGC 1047 

02:40:32.7 

-08:08:51.4 

0.55 

13 

1.8 

1340 

SO-a 

NGC 1052 

02:41:04.8 

-08:15:21.0 

0.71 

32 

1.9 

1510 

E4 

NGC 1060 

02:43:15.0 

32:25:30.0 

0.74 

86 

7.2 

5190 

SO 

NGC 1066 

02:43:49.9 

32:28:29.7 

0.71 

62 

5.8 

4346 

E 

NGC 1361 

03:34:17.7 

-06:15:54.0 

0.68 

50 

6.0 

5272 

E 

UGC 4188 

08:03:24.0 

41:54:53.3 

0.50 

108 

13.9 

9685 

SO 

UGC 4551 

08:44:06.0 

49:47:37.6 

0.31 

13 

2.4 

1749 

SO 

NGC 2675 

08:52:04.9 

53:37:02.3 

0.69 

103 

12.3 

9231 

E 

IC 522 

08:54:34.9 

57:10:00.7 

0.81 

42 

7.1 

5079 

SO 

NGC 2693 

08:56:59.2 

51:20:50.7 

0.69 

76 

6.6 

4942 

E3 

UGC 4702 

08:58:51.2 

38:48:34.4 

0.93 

76 

n.i 

8443 

SO 

NGC 2768 

09:11:37.4 

60:02:14.2 

0.40 

31 

1.7 

1373 

E 

NGC 3265 

10:31:06.7 

28:47:47.7 

0.73 

8 

1.7 

1319 

SO 

NGC 3377 

10:47:42.3 

13:59:09.0 

0.47 

16 

1.0 

665 

E5 

UGC 5928 

10:49:47.1 

51:53:39.0 

0.98 

59 

9.8 

7393 

SO 

NGC 3522 

11:06:40.4 

20:05:07.6 

0.53 

12 

1.1 

1221 

E 

NGC 3539 

11:09:08.8 

28:40:19.9 

0.33 

74 

12.8 

9707 

SO-a 

UGC 6435 

11:25:35.0 

00:46:05.5 

0.78 

99 

10.1 

7604 

SO 

UGC 6683 

11:43:16.2 

19:44:55.4 

0.35 

39 

9.2 

7542 

SO-a 

NGC 3844 

11:44:00.8 

20:01:45.6 

0.26 

46 

8.9 

6779 

SO-a 

NGC 4187 

12:13:29.3 

50:44:29.0 

0.63 

117 

12.1 

9138 

E 

NGC 4458 

12:28:57.6 

13:14:30.9 

0.93 

7 

0.8 

635 

E0 

NGC 4478 

12:30:17.4 

12:19:42.4 

0.81 

15 

1.8 

1349 

E2 

M87 

12:30:49.4 

12:23:28.0 

0.93 

52 

1.8 

1307 

E 

IC 3457 

12:31:51.4 

12:39:25.0 

0.68 

17 

2.1 

1297 

E3 

NGC 4787 

12:54:05.5 

27:04:07.0 

0.32 

52 

10.4 

7585 

SO-a 

NGC 4797 

12:54:55.1 

27:24:45.7 

0.74 

84 

10.1 

7863 

SO 

NGC 4827 

12:56:43.5 

27:10:43.2 

0.79 

61 

9.9 

7630 

E 

NGC 4952 

13:04:58.3 

29:07:19.7 

0.60 

58 

8.2 

5968 

E 

NGC 5004 

13:11:01.5 

29:38:12.1 

0.72 

57 

9.3 

7044 

SO 

NGC 5173 

13:28:25.3 

46:35:29.5 

0.96 

22 

3.2 

2419 

E0 

NGC 5329 

13:52:10.0 

02:19:30.4 

1.00 

102 

9.7 

7109 

E 

UGC 8986 

14:04:15.8 

04:06:44.7 

0.98 

9 

1.8 

1232 

SO 

NGC 5576 

14:21:03.6 

03:16:15.4 

0.66 

30 

2.0 

1487 

E3 

NGC 5638 

14:29:40.3 

03:14:00.5 

0.87 

20 

2.3 

1676 

E 

IC 1024 

14:31:27.1 

03:00:32.2 

0.43 

15 

2.0 

1454 

SO 

IC 1071 

14:54:12.4 

04:45:00.9 

0.63 

107 

10.2 

8302 

so 

NGC 5813 

15:01:11.2 

01:42:07.2 

0.66 

51 

2.6 

1972 

E 

UGC 10261 

16:11:04.0 

52:27:02.1 

0.68 

235 

25.0 

19018 

SO 

Mrk 501 

16:53:52.2 

39:45:36.6 

0.78 

107 

13.3 

10092 

E 

NGC 6364 

17:24:27.3 

29:23:25.4 

0.78 

74 

9.3 

6874 

SO 

UGC 10935 

17:38:43.2 

57:14:21.6 

0.51 

49 

12.7 

8801 

SO 

NGC 7317 

22:35:51.8 

33:56:41.8 

0.95 

56 

8.2 

6599 

E4 

NGC 7432 

22:58:01.9 

13:08:04.4 

0.59 

134 

9.8 

7615 

E 


Notes. 

a Ellipticities at the total major and minor axes are based on the extended 
photometry measurements in the SDSS r-band. 
b Total and half-light radii come from the WISE 3.4 /xm photometry. 
c Radial velocities were obtained from the NASA Extragalactic Database archive 
(http://ned.ipac.caltech.edu/). 

d The Hubble type, according to the Gil de Paz et al. (2007) UV catalog. 


Table 1 includes some of the extended photometry from this 
analysis: ellipticity (r-band) and WISE 3.4 fim total and half- 
light radii. The galaxies span a wide range of total radii from 
7 < R tot < 235 kpc with half-light radii from 0.8-25 kpc. The 
smallest half-light radius (16") is bigger than the point- spread 
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function for the WISE images. The largest galaxies have total 
radii exceeding 100 kpc: UGC 10261 ( 7? tot /7?h = 235/25 kpc), 
NGC 7432 (R t ot/Rh = 134/9.8kpc),andNGC4187(R tot /R h = 
117/12.1 kpc). Approximately 30% of the sample is highly 
elliptical ( a/b < 0.6). The sample has a mean radial velocity, 
cz, of 5396 km s _1 . 

Five of the ETGs in our sample were observed as part of 
the SAURON two-dimensional kinematic survey: NGC 2768, 
NGC 3377, NGC 4458, M87, and NGC 5813. They have been 
classified as slow rotators, indicative of a quiescent stellar 
population (McDermid et al. 2006). Bureau et al. (2011) find 
that the quiescent, slow rotators have blue FUV — NUV colors 
(0.5 < FUV - NUV < 2). We note that the SAURON survey 
analysis is limited to a narrow spatial range within the half-light 
radius. 

3. ANALYSIS 

In this present work, we consider three possible sources for 
the observed UVX in ETGs: (1) a significant population of 
BF1B/EHB stars (f age > 2 Gyr and 7/ if > 14,000 K), where 
EHBs are the combined BHB phase and higher 7/n, (2) a 
substantial population of hot young stars from a recent small- 
scale starburst (£ age < 500 Myr and 1<Z<1.5Z Q , T e g > 
10,000 K), and (3) a significant population of very hot, post- 
main-sequence metal-poor stars (f age > 9 Gyr and Z < 0.5 Z Q ). 
We seek to distinguish between these possibilities by comparing 
their predicted UV/optical/mid-IR colors. 

We use the Flexible Stellar Population Synthesis models 
(FSPS; Conroy et al. 2009; Conroy & Gunn 2010) to create 
composite stellar population (CSP) models. The main advantage 
of the Conroy et al. (2009) models is the inclusion of a BFIB 
phase, and the flexible treatments of the post-AGB and HB 
evolutionary phases. This flexibility allows one to choose the 
age at which the BHB phase is initiated and the percentage 
of thermally pulsating asymptotic giant branch (TP-AGB) stars 
that enter that phase. In essence, adding the BHB /EHB fractions 
to the CSPs drives the colors blueward after £ age = 2 Gyr. 
We combine the FSPS models with the Padova isochrones 
(Marigo & Girardi 2007; Marigo et al. 2008) and the BaSEL 
stellar libraries (Lejeune et al. 1997, 1998), which provide 
mass, age, and metallicity ranges of 0.15 / /M / 100 Mq, 
10 6 ' 6 < t < 10 10 - 2 yr, and 0.005 < Z/Z 0 < 1 .5 (Z Q = 0.019), 
respectively. 

For this analysis, we consider the BHB/EHB phase and 
young star sources of UVX. We also include a comparison 
between subsolar, solar, and super-solar metallicities. We limit 
the metallicities to 0.25, 1, and 1.5 Z Q ; ages between 30 Myr 
and 14.1 Gyr at Alog(t/yr) = 0.025 intervals; apply a Chariot 
& Fall (2000) dust parameter of 0.3; and use the Chabrier (2003) 
initial mass function (IMF). 

We assume that stars on the HB are allowed to become 
BHB/EHB at ages >2 Gyr, based on Yi et al. (1998). They 
show a UV slope increase for 1-1.5 Z Q populations during the 
post-AGB epoch, beginning at ages of >2.5 Gyr. We take three 
different fractions of stars on the HB that are in the BHB (or 
EHB) phase: 0%, 25%, and 50%. Note that, for all populations, 
we have the same total number of stars on the HB ; the fraction 
refers to the division between “ordinary” HB stars and blue (B) 
or extended (E) HB stars. 

The division of stars in the B/E phase between BHBs and 
EHBs is determined by temperature, with higher temperatures 
meaning more E (extended) than B (blue). The increased 
temperature is assumed to be due to small stellar envelope 


masses (<0.05 M Q ), where the radiation is allowed to escape 
freely. This analysis integrates temperature increases of log 
Ay- (log 7/ff) from 0.0, 0.2, and 0.4 dex, which is a measure of the 
change from T e g (e.g., Marigo et al. 2008; Conroy et al. 2009). 
The parameter Ay is the shift in 7/ff for a TP-AGB phase, and 
is included in the FSPS models based on Marigo et al. (2008). 
We chose those values based on the ranges presented for a 
super-solar metallicity isochrone in Figure 1 of Marigo et al. 
(2008) to compare extreme boundaries in our parameter range. 
Finally, we apply three different e-folding times to vary the star 
formation history (SFH) and stellar mass with age (r = 0.2, 0.6 
and 1 Gyr). 

We use the FSPS models to generate tracks for color anal- 
ysis. Specifically, we wish to determine which combination of 
UV/optical/mid-IR color separates populations in as many of 
the following variables as possible: age (f age ), metallicity (Z), 
change in temperature (Ay), SFH (r), and BHB/EHB fraction. 
The plots in Figure 1 provide a sample of the explored color- 
space/parameter combinations. These are an example of three 
color-space combinations with Z = 1 Z Q . The top, middle, and 
bottom rows are NUV -r/g- r, FUV - NUV/NUV - [3.4], 
and FUV — NUV/NUV — r, respectively. The panels in each 
row display the same templates, but the lines are color-coded 
differently to emphasize how the parameters are evolving with 
BHB fraction (left), e-folding time r in Gyr (middle), and the 
increase in effective temperature Ay (logoff)) at HB onset 
(right). For example, the left panels in each row show when the 
BHB fraction takes effect, when the lines have a slower rise in 
NUV — r (top), or when they shift sharply down (blueward) in 
FUV — NUV (middle and bottom); the middle columns high- 
light the track separations for a given r and therefore trace the 
SFH and the right panel does the same for temperature change. 
A more detailed description of the effects of different stellar 
population parameters on the UV, optical, and mid-lR colors is 
included in the Appendix. 

In principle, one can do this with a UV, or UV-optical 
slope in the y-axis, and an optical-near-IR, or optical-mid-IR 
slope in the x-axis. The ultimate goal is to see which color- 
spaces work best, particularly extending the color-spaces into 
the WISE channels (from the near-lR or red optical) to determine 
what combination provides improved discriminatory power. We 
investigated several color-spaces, and determined that a UV 
plus WISE combination does significantly better than the UV 
plus optical (see the Appendix); we address this choice briefly 
in Section 4.2. 

4. RESULTS 

4.1. Radial Color Gradients 

To study radial differences between UV, optical and mid-IR, 
we show color gradients for NUV — r and NUV — [3.4] in 
Figure 2. In Table 2, we list the NUV — [3.4] and NUV — r 
colors at 0.3 and 0.7 R/Rgo fractional radii for each ETG (Rgo 
is the radius at 90% of the total flux in the / -band). For almost 
all galaxies in our sample, there is a trend toward bluer colors 
with increasing radius. To investigate this trend, in Figure 3 we 
plot the histograms of NUV — r and NUV — [3.4] colors at the 
inner (7?; n , R within 50% flux, white) and outer (7? out , R between 
50%-90% flux) radii. The flux-fraction is based on the total 
/--band flux for each galaxy. Both histograms show a statistically 
significant color separation between peaks of approximately 
1 mag. The averages are: NUV — r = 5.7 ± 0.2 and 4.9 ±0.1 
for R[ n and 7? out , respectively, and NUV — [3.4] = 6.1 ± 0.2 
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Figure 1. Flexible Stellar Population Synthesis (FSPS; Conroy et al. 2009), composite stellar population (CSP) templates depicting the evolution of color-space within 
an age range of 0.03-14.1 Gyr. The three rows show NUV — r/g — r (top), FUV — NUV/NUV — [3.4] (middle), and FUV — NUV/NUV — r (bottom) all with 
Z = 1 Zq. The left/middle/right panels in each row show the same templates/color-color values. We color-code the lines differently (left to right, as labeled in the 
top row) to emphasize how the parameters are evolving with BHB fraction, e-folding time r in Gyr, and the increase in effective temperature A t (log(7’ e ff)) at HB 
onset. 


and 5.1 ± 0.2 for R m and R out , respectively. The bluest 
galaxy, Mrk 501, is an outlier in the sample and a known 
blazar. 

We also highlight subpopulations in the histograms based on 
ellipticity and morphological type (Figure 3). Morphological 
type does not appear to be a significant factor in the different 
colors in the inner and outer regions since their colors distribute 
similarly. Highly elliptical systems (b/a < 0.6) have mean R m 
colors that are preferentially reddened by ~2% for NUV — r and 
by ~3% for NUV — [3.4], which is insignificant for our results. 
We address how dust may affect the UV to mid-IR colors in the 
following section. 

4.2. Comparisons to FSPS Templates 

To determine the origins of the color differences, we compare 
the sample UV, optical, and mid-IR colors of the R m and R out 
regions to the CSP models derived from the FSPS templates. In 
Section 3, we described Figure 1 to show how changes in the 
parameters can dramatically alter the paths of the FSPS tracks 
in a given color plot. From this figure and the Appendix, it can 
be seen that the FUV, NUV, r, and [3.4] plots provide a broader 
range in color-space with which to discriminate the dominating 
parameters for a given color region than the other optical color 
combinations. For example, g — r and g — i show a color range 
of approximately a magnitude (Figure 1, and Figures 10 and 11), 


while NUV — r and NUV — [3.4] extend over 6 mag (Figure 1, 
and Figures 12 and 13). Therefore, we focus on these colors in 
the following analysis. 

In Figure 4, we show the FSPS tracks for NUV — r, NUV — 
[3.4], and FUV - NUV versus age (0.2 < t age <14.1 Gyr). The 
left, middle, and right panels for each color-age plot are different 
temperature boosts given to the TP-AGB phase: Ay = 0, 0.2, 
and 0.4 as labeled in the top row. The most significant effects 
on color are metallicity (different colors) and SFH (r). As is 
commonly known, higher metallicities redden the color and 
a smaller r results in a more rapid color increase. A stellar 
population with Z > Z Q and short e-folding time ~0.2 can 
become very red by r age = 2 Gyr. 

The ETG sample colors are shown in gray and orange shaded 
regions that span the R m and R out values, respectively. The darker 
horizontal gray and orange lines mark the mean NUV — r, 
NUV — [3.4], and FUV — NUV values for the corresponding 
color shaded area. As stated previously, the inner radii have a 
significantly redder color than the outer radii in NUV — r and 
NUV — [3.4]. However, this reverses for FUV — NUV, where the 
inner color is, on average, bluer by a small margin (~0. 1 mag). 
As the stellar populations age, the possible physical drivers 
behind the observed colors without assuming priors, quickly 
become degenerate. As a generalized interpretation of the color 
differences in Figure 4, the following results are based on the 
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Figure 2. Color as a function of fractional r-band radii: NUV — r (left) and NUV — [3.4] (right). The fractional radius is taken at R/Rgo, where Rgg is the radius at 
90% of the total flux in the r-band. There is a distinct trend toward bluer colors outward of the half-light radius. The gray shading is the estimated error. Table 2 lists 
NUV - [3.4] and NUV - r values at 0.3 and 0.7 R/Rg 0 - 



NUV— r NUV - [3.4] 

Figure 3. Distributions of NUV — r and NUV — [3.4] colors at the inner ( R m , R within 50% flux, white) and outer (Rout. R between 50%-90% flux, gray) radii. 
The hatches in the top panels accent a subpopulation of highly elliptical galaxies (b/a < 0.6); the bottom highlight SO galaxies. There is a distinct color separation 
between the inner and outer colors, agreeing with the trend seen in Figure 2. This separation occurs at NUV — [3.4] > 6, which is a strong color cut for galaxies with 
low BHB fraction and 1—1.5 Zq (see Section 3 and the Appendix; Figure 13). 


mean colors (gray and orange horizontal lines) for R m and R out 
at ages 2, 5, and 10 Gyr (vertical white lines). 

We summarize the parameter combinations consistent with 
the R in and R mt mean NUV - r, NUV - [3.4], and FUV - NUV 


colors in Table 3. R m tends toward higher metallicities than R out 
for f a g e < 5 Gyr. For both regions, r values are limited to short 
e-folding times (0.2 Gyr) at f age ~ 2 Gyr; the e-folding time 
is degenerate at 5 and 10 Gyr. The BHB fraction is impossible 
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Table 2 

UV-Optical and UV- Mid-IR Colors at Specific Radii 


Object 

NUV 

— 

[3.4] a 

NUV 

— 

[3.4] b 

NUV- 

- r a 

NUV - 

- r ° 


R/Rga ■ 

= 0.3 

R/Rgo - 

= 0.7 

R/Rgo ■ 

= 0.3 

R/Rgo ■ 

= 0.7 

NGC 155 

6.23 

± 

0.03 

5.81 

± 

0.02 

6.09 

± 

0.03 

5.58 

± 

0.02 

NGC 163 

6.37 

± 

0.02 

6.15 

± 

0.02 

6.05 

± 

0.02 

5.77 

± 

0.01 

IC 1698 

6.00 

± 

0.04 

5.14 

± 

0.02 

5.61 

± 

0.04 

4.64 

± 

0.02 

IC 1700 

6.53 

± 

0.03 

6.07 

± 

0.02 

6.17 

± 

0.03 

5.72 

± 

0.02 

UGC 1040 

6.85 

± 

0.05 

6.52 

± 

0.03 

6.76 

± 

0.05 

6.28 

± 

0.03 

NGC 1047 

4.33 

± 

0.01 

4.55 

± 

0.01 

4.38 

± 

0.01 

4.53 

± 

0.01 

NGC 1052 

6.42 

± 

0.01 

6.10 

± 

0.01 

5.89 

± 

0.01 

5.67 

± 

0.01 

NGC 1060 

7.79 

± 

0.03 

7.31 

± 

0.03 

6.88 

± 

0.03 

6.43 

± 

0.03 

NGC 1066 

7.84 

± 

0.06 

7.13 

± 

0.05 

6.91 

± 

0.06 

6.25 

± 

0.04 

NGC 1361 

6.16 

± 

0.03 

5.80 

± 

0.03 

5.93 

± 

0.03 

5.49 

± 

0.03 

UGC 4188 

7.05 

± 

0.05 

6.30 

± 

0.03 

6.75 

± 

0.05 

5.78 

± 

0.03 

UGC 4551 

6.64 

± 

0.02 

6.41 

± 

0.02 

6.53 

it 

0.02 

6.04 

± 

0.01 

NGC 2675 

6.63 

± 

0.04 

6.37 

± 

0.03 

6.64 

± 

0.04 

5.98 

± 

0.02 

IC 522 

6.85 

± 

0.04 

5.93 

± 

0.02 

6.58 

± 

0.04 

5.52 

± 

0.02 

NGC 2693 

6.12 

± 

0.01 

6.05 

± 

0.01 

5.75 

± 

0.01 

5.64 

± 

0.01 

UGC 4702 

6.55 

± 

0.04 

6.05 

± 

0.03 

6.27 

± 

0.04 

5.65 

± 

0.03 

NGC 2768 

6.74 

± 

0.01 

6.26 

± 

0.01 

6.18 

± 

0.01 

5.87 

± 

0.0 

NGC 3265 

3.59 

± 

0.01 

3.66 

± 

0.01 

3.34 

± 

0.01 

3.39 

± 

0.01 

NGC 3377 

6.18 

± 

0.01 

5.57 

± 

0.01 

5.75 

± 

0.01 

5.33 

± 

0.0 

UGC 5928 

6.09 

± 

0.05 

5.78 

± 

0.03 

6.09 

± 

0.05 

5.58 

± 

0.03 

NGC 3522 

5.60 

± 

0.02 

5.12 

± 

0.01 

5.90 

± 

0.02 

5.17 

± 

0.01 

NGC 3539 

7.36 

± 

0.09 

6.49 

± 

0.04 

7.00 

± 

0.09 

5.96 

± 

0.04 

UGC 6435 

6.46 

± 

0.03 

5.96 

± 

0.02 

6.31 

± 

0.03 

5.63 

± 

0.02 

UGC 6683 

6.59 

± 

0.05 

6.42 

± 

0.04 

6.61 

± 

0.05 

6.17 

± 

0.03 

NGC 3844 

6.91 

± 

0.05 

6.54 

± 

0.03 

6.72 

± 

0.05 

6.20 

± 

0.03 

NGC 4187 

6.81 

± 

0.03 

6.35 

± 

0.02 

6.25 

± 

0.02 

5.93 

± 

0.02 

NGC 4458 

5.88 

± 

0.01 

5.47 

± 

0.01 

5.77 

± 

0.01 

5.41 

± 

0.01 

NGC 4478 

6.38 

± 

0.01 

6.07 

± 

0.01 

6.17 

it 

0.01 

5.85 

± 

0.01 

M87 

5.55 

± 

0.01 

5.65 

± 

0.01 

5.19 

± 

0.01 

5.35 

± 

0.01 

IC 3457 

4.62 

± 

0.01 

4.56 

± 

0.01 

4.85 

± 

0.01 

4.71 

± 

0.01 

NGC 4787 

6.56 

± 

0.05 

5.97 

± 

0.03 

6.53 

± 

0.05 

5.70 

± 

0.03 

NGC 4797 

6.61 

± 

0.01 

6.25 

± 

0.01 

6.04 

± 

0.01 

5.79 

± 

0.01 

NGC 4827 

6.20 

± 

0.01 

6.04 

± 

0.01 

6.05 

it 

0.01 

5.68 

± 

0.01 

NGC 4952 

6.28 

± 

0.02 

6.11 

± 

0.02 

6.19 

± 

0.02 

5.75 

± 

0.01 

NGC 5004 

6.37 

± 

0.02 

6.06 

± 

0.01 

6.24 

± 

0.02 

5.69 

± 

0.01 

NGC 5173 

4.61 

± 

0.01 

4.31 

± 

0.01 

4.67 

± 

0.01 

4.23 

± 

0.01 

NGC 5329 

6.29 

± 

0.01 

5.99 

± 

0.01 

6.01 

± 

0.01 

5.68 

± 

0.01 

UGC 8986 

5.19 

± 

0.01 

4.81 

± 

0.01 

5.32 

± 

0.01 

4.96 

± 

0.01 

NGC 5576 

6.05 

± 

0.01 

5.71 

± 

0.01 

5.86 

± 

0.01 

5.52 

± 

0.01 

NGC 5638 

6.35 

± 

0.01 

5.96 

± 

0.01 

5.99 

± 

0.01 

5.71 

± 

0.01 

IC 1024 

4.75 

± 

0.02 

4.32 

± 

0.01 

3.64 

± 

0.01 

3.64 

± 

0.01 

IC 1071 

6.98 

± 

0.02 

6.51 

± 

0.02 

6.37 

± 

0.02 

5.98 

± 

0.02 

NGC 5813 

6.47 

± 

0.01 

6.14 

± 

0.01 

6.00 

± 

0.01 

5.72 

± 

0.01 

UGC 10261 

6.63 

± 

0.02 

6.21 

± 

0.01 

6.09 

± 

0.01 

5.65 

± 

0.01 

Mrk 501 

2.56 

± 

0.01 

3.09 

± 

0.01 

2.18 

± 

0.01 

2.46 

± 

0.01 

NGC 6364 

6.71 

± 

0.02 

6.27 

± 

0.01 

6.47 

± 

0.02 

5.91 

± 

0.01 

UGC 10935 

6.92 

± 

0.08 

5.89 

± 

0.04 

6.39 

± 

0.08 

5.43 

± 

0.04 

NGC 7317 

6.96 

± 

0.02 

6.78 

± 

0.02 

6.74 

± 

0.02 

6.18 

± 

0.02 

NGC 7432 

6.63 

± 

0.03 

6.29 

± 

0.02 

6.21 

± 

0.03 

5.82 

± 

0.02 


Notes. 

a Colors for NUV — [3.4] and NUV — r measured at the r-band fractional radii 
R / 7? 90 = 0.3. See Figure 2. 

b Colors for NUV — [3.4] and NUV — r measured at the r-band fractional radii 
R/Rgo = 0.7. See Figure 2. 

to isolate at any temperature or age in these plots. However, 
the inner mean color is consistent with having little to no BHB 
population. 

Additionally, it appears that temperature is only important 
at the extreme end A j = 0.4. Despite the degeneracies, we 
find indications that age is primarily responsible for the color 
difference in our ETG sample. From just the three sampled 


ages, we show the progression of the different modeled stellar 
populations that might lead to color differences. 

To explore the premise that age is the main cause of the color 
differences in the ETGs, we take isochrones of the synthetic 
photometry in Figure 4 at f age = 2, 5, and 10 (white vertical 
lines in Figure 4), and plot them in FUV — NUV, NUV — r, 
NUV — [3.4] color-space. We do this in Figure 5 by plotting 
FSPS lines on top of the gray-scale contours that include both 
the R m and R out colors of the ETG sample. We outline the 
distribution of colors for R out with the dashed contour lines. 
We also add dust vectors (orange arrows crossing the Z = 1 
and 1.5 Z Q tracks) with an increasing dust parameter from 0-1 
using the Chariot & Fall (2000) dust models. The top, middle, 
and bottom rows are the isochrones at ages 2, 5, and 10 Gyr as 
labeled. 

Our first assessment is that the density peaks are clearly 
separated, and dust may be ruled out simply by observing the 
direction of the vectors in Figure 5. If dust was preferentially 
reddening the inner half-light of the ETGs, the FUV — NUV 
R[ n colors would lie redward. We do not mean that these are 
entirely dustless (we actually include dust in the FSPS models; 
see Section 3). Rather, we assume that dust is not a strong 
enough factor to cause the color difference. 

The 2 Gyr plots (top row in Figure 5) indicate that the colors 
for the inner regions could only be explained by super-solar 
metallicities >1.5 Z Q . At the same age, the outer regions are 
more likely tohosta 1-1.5 Z Q stellar population with a moderate 
BHB fraction. At 5 Gyr (middle plots), the R n] colors could 
be caused by higher metallicity/moderate BHB fractions and 
the outer regions would have a lower metallicity and higher 
numbers of BHBs. At 10 Gyr, the isochrones give narrower 
ranges and the metallicity dependence is stark: A’ ]n would be 
dominated by a 1-1.5 Z Q stellar population and an insignificant 
BHB fraction; R out would be dominated by a < 1 Z Q stellar 
population and may include slightly more BHBs. We also note 
that for 2-5 Gyr the dominating SFH is r = 0.2-0. 6; at 10 Gyr 
the range is widened to all r. 

The above conditions indicate two possible scenarios: (1) if 
the bulk of the bulge and disk mass coevolved, the metallicities 
must be different by ~1 Z Q and (2) the ETGs formed in an 
inside-out process with at least two major stages of growth 
>1 Gyr apart. These appear to be two distinctly different 
mechanisms of formation, and we deepen our analysis in the 
following section to determine the likelihood of either scenario. 

4.3. Age Estimates for R, n and R out 

To constrain the likelihood of the results based on Figure 5, 
we estimate the ages for the inner and outer regions of the 
ETGs. However, we showed in the previous section that the 
high degeneracy of the synthetic photometry makes it difficult 
to derive accurate ages for each of the objects (e.g., Figure 4). In 
order to check the results, given the degeneracies, we resample 
the data to find the most likely statistical fit to the ages given a 
fixed set of parameters. We fit the data to the FSPS templates, 
with the FUV - NUV, NUV - r, and NUV - [3.4] colors. We 
resampled the inner and outer colors by replacement estimates 
(bootstrapping) within the errors for 500 iterations to estimate 
errors and probabilities for the y 2 fitting. The y 2 fitting was 
done between the resampled observed colors and the FSPS 
data. After minimizing the y 1 for each iteration, we used 
the maximum likelihood for each parameter combination. The 
probability distribution functions were used to estimate the 
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tage (Gyr) tage (Gyr) t a ge (Gyr) 

Figure 4. To illustrate the significant degeneracies with age, we plot color vs. age from 0.2 < t a ge < 14.1 Gyr, with the ETG sample (orange and gray shading), and 
the FSPS colors (NUV — r, NUV — [3.4], FUV — NUV) with temperature boosts A j = 0, 0.2, and 0.4 from left to right (as labeled). The magenta, black, and blue 
lines are Z = 0.25, 1.0, and 1.5 Z©, respectively. The solid, dashed, and dotted lines represent the 0, 0.2, and 0.5 BHB fractions initiated at t a ge = 2 Gyr. Each line 
style has three lines (not labeled) for the three star formation histories (r = 0.2, 0.6, 1.0). The range of color for the inner and outer regions of the ETGs are shaded 
in gray and orange bands, respectively. The solid horizontal lines indicate the mean values for the respective regions. The white vertical lines are placed at 2, 5, and 
10 Gyr, which we inspect in Figure 5. 


Table 3 

Possible Parameter Solutions for the Mean R m and R ou t Colors at 2, 5, and 10 Gyr 


fage 


A r =0 



Aj- = 0.2 



o 

II 

<1 


Z(Z 0 ) 

r (Gyr) 

/bhb 

Z(Z 0 ) 

r (Gyr) 

/bhb 

Z(Z 0 ) 

r (Gyr) 

/bhb 

^in 

2 Gyr 

1.5 

0.2 

<0.25 

1.5 

0.2 

<0.25 




5 Gyr 

1-1.5 

0.6 

0-0.5 

1-1.5 

0.2-0.6 

<0.25 

1-1.5 

0.2-0.6 

<0.25 

10 Gyr 

0.25-1.5 

0.2-1 

0-0.5 

1-1.5 

0.2-1 

<0.25 

1-1.5 

0.2-1 

<0.25 

^OUt 

2 Gyr 

1-1.5 

0.2 

0-0.5 

1-1.5 

0.2 

>0.25 

1-1.5 

0.2 

0-0.5 

5 Gyr 

0.25-1.5 

0.2-0.6 

0-0.5 

0.25-1 

0. 2-0.6 

0-0.5 

0.25-1 

0.2-0.6 

0-0.5 

10 Gyr 

<1 

0.2-1 

0-0.5 

0.25-1 

0.2-1 

0-0.5 

0.25-1 

0.2-1 

0-0.5 


weighted average age. The ages and errors were calculated for 
each BHB, r, Z, and A T combination for each object. 

In Figure 6 we show the effects that BHB fraction, metallicity, 
and SFH have on the estimated ages (f age (in) and f age (out)). Each 
data point is a weighted average of the entire ETG sample, 
testing different parameter configurations (e.g., one point for 
/ BHB = 0, r = 0.2, and Z = 1 Z Q ). We do not separate the 
temperature boosts here since Figure 4 indicates that the effects 
on color are negligible for this analysis. 

The points reflect the assumption that both the inner and 
outer regions are evolving with the same set of parameters. We 
are essentially testing one condition for the first scenario in 
Section 4.2. For example, the magenta square (Z = 1.5 Z Q , 
/bhb = 0.5, and r = 1), easily identified as the outlier, is 


fixing the prior that the parameter combinations are the same 
for the inner and outer regions. In other words, each point is the 
t age coordinate for the inner and outer regions having the same 
metallicity, BHB fraction, and r. With this assumption, we find 
that nearly all points are below the equality line, and many 
(12 points) have Af age > 1 Gyr. 

In order for the inner and outer regions to have equal ages, 
they must have a certain set of properties. Over half of the 
points within errors of the equal age line have low metallicity 
(double circles for Z = 0.25 Z Q ) and some fraction of 
BHBs (diamonds and squares). This constrains our hypothetical 
scenarios mentioned in the previous section. 

Average ages derived from low metallicity (0.25 Z Q ) show 
the least difference between f a ge(i n ) and t age (out). The average 
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NUV-r NUV — [3.4] 


Figure 5. FUV — NUV, NUV — r, and NUV — [3.4] plots with the distribution of inner/outer colors, and isochrones at 2 (top), 5 (middle), and 10 (bottom) Gyr. 
The gray-filled contours include the inner and outer radii colors, while the outer regions are delineated with dashed contour lines. The peaks of the two regions are 
clearly separated in color-space. The dashed, dotted, and solid FSPS lines are Z = 0.25, 1, and 1.5 Zq, respectively. The temperature increases are At = 0 (black), 
0.2 (magenta), 0.4 (cyan). All lines include SFH ^-folding times r = 0.2, 0.6, and 1 Gyr. The BHB fraction at points along the isochrone tracks are marked by the 
symbols as labeled. We include dust vectors (orange) with an increasing dust parameter 0-1, using a Chariot & Fall (2000) dust model. 


age separation is 0.28, 1.13, and 0.9 Gyr for Z = 0.25, 1, 
and 1.5 Z Q , respectively. When selected by BHB fraction, the 
average differences are 1.21, 0.86, and 0.24 Gyr for /bhb = 0, 
0.25, and 0.5, respectively; r shows little difference 0.83, 0.8, 
and 0.69 Gyr for r = 0.2, 0.6 and 1, respectively. 

Given these priors, 67% of the ages for the inner regions are 
older (to the right of the shaded region), and have an average 
age difference of 1 . 1 Gyr. The mean difference for all points is 
0.8 Gyr with a range of —0.3 to 1.9 Gyr if the parameters are 
assumed to be equal. The average age for the inner region is 
7.0 ± 0.3 Gyr and the outer is 6.2 ± 0.2 Gyr. The assumptions 
we used to derive these results indicate that the outer regions 
are likely to have formed 0.5-1. 2 Gyr after the inner regions. 

4.4. Age-Metallicity Ratios 

We explore the possibility of an age-metallicity difference 
between R [n and R out in our sample. Since we are limited to 
three metallicities, the results we present are a crude estimate 


of how metallicity and age behave between these regions. 
Figure 7 shows the inner to outer average age ratios for different 
■Zirmer and Z outer metallicity combinations. The error bars span 
the range of ages for the set of parameters at that metallicity 
ratio (e.g., BHB, t). We fit the Alog(Age) — Alog(Z) values 
with the following result: 

Alog(Age) = — 0.16Alog(Z). 

Figure 7 shows that, for Zj nner /Z outer < 1, it is more likely 
that the stellar populations in the inner half-lights of our sample 
are older than the outer even though the range of ages spans 
widely about the mean. For Zi nner /Z outer > 1, the ages span a 
narrow range around tinner/ Guter = 1. For metallicity gradients 
of decreasing Z with radius, the estimated age difference for 
the ETGs increases, resulting in an average of 1.4 Gyr (older 
toward the centers), which is nearly double the average when 
we assumed Zinner/Zouter = 1 (see Section 4.3). This further 
supports the hypothesis that the color differences are likely due 
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dinner (Gyr) 

Figure 6. Estimated ages in Gyr for R m and R out considering the different CSP 
parameter combinations, assuming homogeneity of the stellar populations in the 
inner and outer regions. The ages are combined weighted averages for the inner 
and outer regions for the 49 ETGs in our sample (one symbol is the average for 
all ETGs). The different triangles, diamonds, and squares represent /bhb = 0, 
0.25, and 0.5, respectively. The black, blue, and magenta colors are r = 0.2, 
0.6, and 1, respectively. Symbols on double circles and double diamonds have 
Z = 0.25 and 1 Zq, respectively, while all other symbols have Z = 1.5 Zq. 
The dotted lines trace metallicities along a fixed /bhb and r with metallicity. 
The shaded region highlights the equal age line within the average error. The 
average error is illustrated by the error bars in the upper left comer. We include 
equal age ±1 Gyr dashed lines for gauging the age differences. 



Z inner / Z 0 uter 


Figure 7. We plot the estimated age ratios for each metallicity combination 
between the inner and outer regions of the ETGs. The solid line is the fit of the 
ETG sample (diamonds). The error bars show the range of ages for the different 
BHB, r combinations. The horizontal gray area shades the equal age line within 
a 10% margin. 


to multi-stage evolution and less of a metallicity difference. 
However, the information presented here is not conclusive of a 
presence, or lack thereof, of a metallicity gradient. 

4.5. BHBs and the Color NUV — [3.4] Color Difference 

The majority of our analysis is dedicated to determining 
whether the color difference is indicative of co-evolving or 
multi-stage evolution. We also want to know the source of UV 
emission in these ETGs since we have shown our sample is 
unlikely to have undergone recent star formation. The age esti- 
mates indicate that significant UV emission is largely produced 



Figure 8. We show the stacked distributions of A(NUV — [3.4]) determined by 
imposing different BHB fractions at the inner and outer radii as labeled. For 
example, the dark gray labeled 0/0.25 refers to the color at /bhb = 0 and 0.25 
for R[ n and R out , respectively. The models with higher BHB fraction in the outer 
radii are more likely to have a 1 mag color difference (light blue, dark gray, and 
light gray). 


by sources in an older stellar population (f age > 2 Gyr; see 
Figures 5 and 6). One question we asked is whether a sig- 
nificant fraction of BHB stars can accomplish this. Can the 
overall trend of bluer NUV — [3.4] and NUV — r with in- 
creasing galactocentric radius be due to the presence of >25% 
BHBs in a 2-5 Gyr old stellar population? For example, in 
Figure 5, at 10 Gyr the BHB fraction is important only within 
a small area of color-space: 0.6 < FUV — NUV < 1.5, 
4.4 < NUV - r < 5.7, 4.2 < NUV - [3.4] < 6.3. The mean 
R in /R out values are: (FUV - NUV) = 1.2/1. 3, (NUV - r) = 
5. 7/4. 9, and (NUV- [3.4]) = 6.1/5. 1, which do fall into these 
ranges of color-space, indicating that BHBs could be a source 
of the UV emission. 

To investigate how they may contribute to the 1 mag color 
difference between R m and R out , we assign different BHB 
fractions between the inner and outer radii and extract the 
corresponding NUV and [3.4] photometry from FSPS templates 
using the estimated ages for the ETG sample at R; n and R out 
(see Section 4.4). Figure 8 shows the stacked distributions of 
A(NUV - [3.4]) = (NUV - [3.4]) inner - (NUV - [3.4]) outer ). 
Higher /bhb in the outer regions are most likely to make colors 
redder in the centers. In this case, a jump of /bhb = 0.25 in the 
outer region produces this effect. 

4.6. Mass, Star Formation, and UV-Mid-IR Color 

We explore whether the colors in the inner and outer radii are 
related to star formation rates and stellar mass (SFR and M*, 
respectively) derived from integrated photometry (e.g., for low- 
surface-brightness or high-redshift sources). We first compute 
SFRs and stellar masses from the integrated magnitudes (for 
the whole system in each case) using the UV calibration of 
Kennicutt (1998, Equation (1)) and the 2MASS K s versus g — r 
calibration of Bell et al. (2003). 

In Figure 9 we plot these values against the NUV — [3.4] 
colors for the individual sources and find that (1) the SFRs 
show a scatter for both the inner and outer colors and (2) there 
is an increase in stellar mass with redder color. This increase in 
mass with more red ETGs is consistent with other observations 
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SFR (M 0 /yr) LOG (NT /M 0 ) 

Figure 9. NUV — [3.4] inner (gray pentagons) and outer (stars) region colors as a function of star formation rate (left) and stellar mass (right). The black circles and 
dotted line show the binned averages of the total flux NUV — [3.4] colors for SFR and M* (ASFR = 0.05 Mq yr -1 and Alog (M* /Mq) = 0.2). As expected, the ETGs 
are quiescent, and there appears to be no correlation between the UV-mid-IR color and SFR. The average error for NUV — [3.4] is 0.07 mag. The symbols are larger 
than the photometric errors. 


(e.g., Bell et al. 2003). The binned averages are the total flux 
colors (not inner and outer but the entire galaxy), averaged 
over SFR and stellar mass bins. The SFR is very low for all 
of the galaxies; therefore, we can assume the bluer colors for 
the outer radii have not been created by recent star formation. 
The majority of the ETGs are within the massive galaxy range 
(>10 n Mq). We note that the Bell et al. (2003) analysis does 
not include contributions of HB stars to the UV. Flowever, it is 
useful to discern if the color difference could be caused by a 
primarily young stellar population. 

5. DISCUSSION 

We have shown that a clear difference in the UV to IR color 
between the inner and outer regions of our sample of 49 E/S0 
galaxies exists. Moreover, we see no significant variation be- 
tween the E and SO Hubble types. In this section, we discuss 
which stellar mass assembly histories are consistent with these 
observations based on our analysis in the previous section. The 
CSP predictions of color based on BHBs, ages, and metallicities 
illustrate that secular coevolution of a bulge and disk progen- 
itor would be difficult to explain without a significant internal 
mixing mechanism. 

The first result we found was the 1 mag difference between 
the colors measured for the inner half-light and the outer 
radii. We hypothesized that a different stellar population must 
be dominating the respective regions since color gradients in 
ETGs and globular clusters have been repeatedly measured. For 
example, Carter et al. (201 1) use GALEX and 2MASS to discuss 
the radial UVX in ETGs; they find a steep gradient between the 
inner and outer regions with an increase in UVX at the core. 
They attribute this to increased a enhancement from Type II 
SNe, and perhaps He abundance, rather than to dry mergers 
in an outside-in scenario. Moreover, recent papers by Forbes 
et al. (201 1) and Arnold et al. (201 1) present optical color radial 
profiles for a nearby globular cluster (NGC 1407 and SO galaxy 
NGC 3115, respectively). Both show evidence of a two-stage 
formation history using color analysis with g — i. They find a 
color difference and decrease in metallicity [Fe/H] between the 
inner and outer regions. 


We took advantage of the better-defined separations of the 
FSPS tracks in the NUV — [3.4] color-space, and augmented 
this with the NUV — r and FUV — NUV color to narrow 
the possibilities between the well-known degeneracies. In the 
density plots (Figure 5), the separate peaks for inner and outer 
fall in regions depending on the assumed age, and support 
the idea that EHBs are the possible source of UV emission. 
Furthermore, this is consistent with being driven primarily by 
differences in age with secondary contributions from differences 
in BHB fraction, metallicity, and SFH. 

Based on these assessments, and as we mentioned in 
Section 4.2, two possible formation scenarios are consistent 
with comparison of the FSPS colors and the different colors of 
the inner and outer regions in the ETG sample: (a) if the bulge 
and disk coevolved, the metallicities must be different by > 1 Z G 
and (b) the ETGs formed in an inside-out process with at least 
two major stages of growth >1 Gyr apart. 

For the first case, the inner parts of the galaxy would need to 
exhibit a strong metallicity gradient within a short time-frame 
through self-enrichment or merging. If the hierarchical model 
accurately describes the formation of the bulk of the bulge mass, 
then the timescales for relaxation are on the order of 0.1-1 Gyr 
(e.g., Barnes & Hernquist 1996), and pseudobulge formation at 
timescales >2 Gyr. Surveys of globular clusters in M31 show a 
mix of metal-poor to metal-rich clusters across the galaxy with 
metal-rich clusters predominately in the central 10 kpc radius 
(Saito & lye 2000; Perrett et al. 2002). Saito & lye (2000) model 
the timescale to produce a suitable metallicity gradient over the 
galactocentric radius, and calculate that self-enrichment with 
a collapsing disk would take ~5 Gyr, while a merging event 
would take ~2 Gyr. Both surveys conclude that merging with 
low-metallicity galaxies is most likely to have produced the 
metallicity gradients, supporting multi-stage evolution. If this is 
typical for most ETGs, evidence of a second-stage or multi-stage 
event would be expected. The timescales of self-enrichment are 
too long for the inner and outer stellar populations to have 
coevolved, which leads us to the second scenario. 

Our results are more consistent with an inside-out multi- 
stage evolution. Whether or not the color difference between 
the inner and outer radii is caused by a metallicity gradient, the 
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centers of the ETGs must be on average ~ 1 Gyr older, according 
to our results. The question remains what merger process, 
gas-rich or gas-poor, is most likely to produce the age/color 
differences observed. 

In principle, several processes can contribute to stellar mass 
assembly in local massive galaxies. First, in situ star formation 
may be triggered by accretion of gas from the intergalactic 
medium (IGM). Mergers may also play a part through “wet” 
(dissipational) mergers between gas-rich galaxies, or “dry” 
(dissipationless) mergers between gas-poor galaxies. The stellar 
ages for wet mergers tend to be younger, supplying the gas 
for starbursts. A “frosting” of gas-rich mergers is often used 
to explain the dispersion in ages within ETGs (Trager et al. 
2000; Sanchez-Blazquez et al. 2006, 2009). For dry mergers, 
the most massive bulge-dominated ETGs form by the accretion 
of galaxies with older stellar populations onto the central bulge. 

Each of these may occur via major or minor mergers, and 
observational evidence varies on which is most influential in 
terms of mass building (e.g., van Dokkum 2005; Cox et al. 2006; 
Sanchez-Blazquez et al. 2009). The analysis in van Dokkum 
(2005) shows that 71% of bulge-dominated galaxies have some 
signature of interaction (i.e., tidal disturbances), and 35% have 
had a recent dry major merger. Sanchez-Blazquez et al. (2009) 
did a spectral analysis of the van Dokkum (2005) sample and find 
that the sample easily divides into a range of wet to dry merger 
stages in age and metallicity. The metallicities are 1-1.5 Z Q for 
all merger scenarios (Sanchez-Blazquez et al. 2006, 2009), but 
the dry mergers tend toward slightly higher metallicities and 
ages >5 Gyr. Simulations by Cox et al. (2006) show that wet 
mergers leave remnants that are smaller, with higher rotation, 
and a more disk-like shape. Boxy, massive ellipticals can result 
from dry mergers, but the observed nearby ETGs are likely 
formed by a combination of wet/dry merging. Our sample 
consists of more massive ellipticals, and the estimated ages 
(~7 Gyr) support dry merging for z < 0.3. If a metallicity 
gradient exists, dry merging with higher metallicity galaxies 
may be consistent with our age-Z trends. We find inconsistent 
trends with our first scenario, where we proposed that higher 
metallicities at the center would cause the red colors. However, 
our calculations will need to be verified with a more rigorous 
age-Z radial profile fitting. 

There is evidence that the contribution from wet and dry 
mergers may result in a distinct two-stage formation history 
for ETGs. Observations show that the inner and outer regions 
of nearby galaxies decouple at some point, implying that a wet 
process is producing transitions and velocity dispersions indica- 
tive of recent star formation (Treu et al. 2005; Emsellem et al. 
2007; Kormendy et al. 2009; Arnold et al. 2011; Forbes et al. 
201 1). In an inside-out evolution, the oldest stellar populations 
exist in the center, and the outer halo or disk has undergone 
more recent star formation via minor gas-rich mergers. 

Recent observations from Nelson et al. (2012) support inside- 
out evolution for rapidly forming disks at z ~ 1, and agree 
with our age and multi-stage formation time-frames. Their 
results show an extended Ha disk, outside of an older central 
bulge, compared with the A 1 - band stellar continuum (r,,(Ha) > 
1.3 r e (R)). Damjanov et al. (2009) compare the restframe optical 
radii of compact passive galaxies at 1 < z < 2 with nearby 
ETGs. They show that the inner cores of the nearby ETGs are 
consistent with the sizes of the higher redshift compact galaxies. 

Our results are compatible with a two- or multi-stage, inside- 
out formation history, in which the core formed through rapid 
star formation via wet mergers or cold gas accretion from the 


IGM at least ~7 Gyr ago (z >2). After ~l-2 Gyr (at z ~ 1), 
a second phase of minor-merging would induce a new major 
starburst phase possibly in an outer disk of the galaxy. Our age 
estimates require that the inner and outer regions are both greater 
than 2 Gyr old, and differ by an average of 1 .4 Gyr. 

Additionally, the outer regions are consistent with having an 
elevated fraction of EHBs that can begin as early as 2-5 Gyr, 
based on the synthetic photometry (see Figure 4). Our results 
show that an increase of 0.25 BHB fraction in the outer regions 
can produce the 1 mag color difference we observe. The 
BHB/EHB phases are predicted to have lifetimes of ~10 Myr, 
so a population of stars in this phase would be very short- 
lived compared to the average age of >4 Gyr in the inner 
regions. 

A possible alternative explanation is differential extinction 
between inner and outer regions. However, the extinction vectors 
(orange arrows) shown in Figure 5 imply they would need 
to be over a magnitude of extinction, and the FUV — NUV 
would show the opposite separation between the inner and outer 
regions, which we do not observe. It is possible but unlikely, 
though we cannot completely rule it out from our data. This 
may explain why we see a slight bias for the highly elliptical 
galaxies ( b/a < 0.6) at the inner radii (see Figure 3). 

6. CONCLUSIONS 

The longstanding question of what causes UVX is likely not 
answered by one simple population (old or young) or parameter 
tweak. We explore the origin of the UV emission in our sample 
using stellar population synthesis models. We also chose within 
a range of values, the parameters (BHB fraction, Z, t, At) that 
best fit the UV, optical, and mid-IR colors of our sample. We 
then looked at the spatial breakdown of UV to optical and UV 
to mid-IR colors. Through this approach, the following results 
were observed. 

1. WISE and GALEX colors FUV-NUV and NUV- [3.4] are 
highly effective at separating the parameters that drive the 
observed colors. Our CSP tracks show that, at 10 Gyr, the 
BHB fraction is important only within a small area of color- 
space: 0.6 < FUV - NUV < 1.5, 4.4 < NUV - r < 5.7, 
4.2 < NUV - [3.4] < 6.3. 

2. The 49 ETGs in this sample exhibit a strong color difference 
with bluer colors on the outside of the galaxies in both 
UV-optical and UV-mid-IR. We extract the photometry 
for the inner half-light and outer 50%-90% and find a clear 
color difference independent of E/S0 type. The average 
value of the ETG R ia /R 0 u\. colors are NUV — r = 5. 7/4. 9, 
and NUV - [3.4] = 6.1/5. 1. 

3. Increasing the fraction of BHBs in the outer regions by 0.25 
can create the observed color difference, where the inner 
are redder than the outer regions by 1 mag in NUV — [3.4]. 

4. We find that the properties of the inner and outer regions 
are significantly different, either in age, metallicity, and/or 
the existence of BHBs. We discuss two formation scenarios 
based on our results: (1) if the bulge and disk coevolved, 
the metallicities must be significantly different (0.25 Z Q in 
the outer regions, > 1 Z Q in the centers) and (2) the ETGs 
formed in an inside-out process with at least two major 
stages of growth > 1 Gyr apart. Our age estimates indicate 
that the second scenario is most likely. 

5. The average ages are estimated to be 7.0 ± 0.3 Gyr (inner) 
and 6.2 ± 0.2 Gyr (outer), with a minimum of 2.6 Gyr. Even 
when we assume homogeneity of parameters over radius, 
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Figure 10. Flexible Stellar Population Synthesis (FSPS; Conroy et al. 2009), composite stellar population (CSP) templates depicting the evolution of NUV — r vs. g — 
r color-space within an age range of 0.03-14.1 Gyr. The three rows depict tracks for different metallicities Z = 0.25, 1.0, and 1.5 Z© as labeled. Each row/metallicity 
shows the same templates/color-color values. The lines are colored differently (left to right) to emphasize how the parameters are evolving as labeled in the top row 
panels: BHB fraction, e-folding time r in Gyr, and the increase in effective temperature A j (log(r e ff)) at HB onset. 


there is evidence of multi-stage evolution where the outer 
regions are likely to have formed at least ~0.8 Gyr after the 
inner regions (—0.3 < (A f age } < 1.9 Gyr). 

6. The age differences estimated for metallicity gradients with 
an increase in Z at larger radial distances are younger at the 
center by an average of ~300 Myr. For increasing Z with 
radius they present age differences (t age (in)— t age (out)) of 
1.4 Gyr with a range of 0.6-2. 2 Gyr. This suggests that age 
and other stellar population properties are contributing to 
the color difference rather than metallicity. 

7. Since the estimated ages are beyond the lifetimes of star- 
forming regions to contribute significantly to the UV 
emission, we assume that BHBs or EFIBs are the primary 
source of the UVX, agreeing with previous results (e.g., 
Greggio & Renzini 1990). The average colors fall within 
the ranges predicted for BHB fractions greater than 0.25 in 
NUV — [3.4] in the outer regions. 


These combined results lead us to the conclusion that the 
UV observed in these massive ETGs is likely caused by an 
extreme HB phase of older stars and that the ETGs went through 
a multi-stage evolution that is coincident with an inside-out 
cessation. We do not discount that star formation is causing UV 
in other ETGs, but our particular selection may have biased 
us toward more massive ETGs that are slow rotators and not 
star-forming. A larger sample will be explored in detail in a 
subsequent paper to look at the individual properties of the 
galaxies and their environments in more detail. Previous studies 
do not probe the range of ETG type with the radially binned 
extended photometry, and these very extended objects would 
not have the coverage out to large radii that WISE naturally 
provides. 
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Figure 11. FSPS CSP templates depicting the evolution of NUV — g vs. g — i color-space from ages 0.03-14.1 Gyr. See Figure 10 for a full description. 
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APPENDIX 

COLOR ANALYSIS OF UVX IN STELLAR POPULATIONS 

We used the FSPS photometry to determine the best approach 
for comparing U V-optical, optical-optical, optical-mid-IR, and 
UV-mid-IR colors. The analysis was carried out as described in 
Section 3. We briefly recap the CSP model parameters here. We 
combine the FSPS models with the Padova isochrones (Marigo 
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Figure 12. FSPS CSP templates depicting the evolution of FUV — NUV vs. NUV — r color-space from ages 0.03-14.1 Gyr. See Figure 10 for a full description. 


& Girardi 2007 ; Marigo et al. 2008) and limit the metallicities 
to 0.25, 1, and 1.5 Z G , ages between 30 Myr and 14.1 Gyr at 
A Iog(f/yr) = 0.025 intervals, apply a Chariot & Fall (2000) 
dust parameter of 0.3, and use the Chabrier (2003) IMF. 

For BHBs, we assume that stars on the HB are allowed 
to become BHB/EHB at ages >2 Gyr, taking three different 
fractions of stars on the HB that are in the BHB (or EHB) 
phase: 0%, 25%, and 50%. To address the EHB phase, assuming 
these are BHBs with higher effective temperatures, we include 
temperature boosts Ay = 0.2, 0.4 dex. 

In Figures 10-13, we show how colors evolve with changes 
in the parameters, using NUV — r /g — r, NUV — g/g — i , 
FUV - NUV /NUV - r , FUV - NUV /NUV -[3.4], Each panel 
from left to right shows the same FSPS tracks, but the lines are 
colored differently based on parameter values to emphasize the 
dependence of that parameter on the position in color-space (see 
labels in the top row). For example, the top row in Figure 10 
shows (left) the BHB fractions of 0, 0.25, and 0.5 depicted in 
black, blue, and brown, respectively; (middle) the SFH r = 0.2, 


0.6, and 1 Gyr in black, blue, and brown, respectively; (right) 
the change in effective temperature Ay at 0, 0.2, and 0.4 dex 
in black, blue, and brown, respectively. The consecutive rows 
show how the FSPS tracks change with Z as indicated in each 
panel. 

Separation in color-space for each variable is not strong for 
g — r, NUV — g, and g — i. All three have smaller ranges 
and smaller separations, leading to significant degeneracies. 
Concerning the former point, the color ranges in the panels 
in Figure 10 are very narrow, spanning 0.7 in optical, compared 
with the NUV — r and NUV — g (~6). For the second point, 
the largest difference between the BHB and r lines for g — i 
and g — r are <0.1, and 0.3, respectively. Separation in r is 
reasonable, but separation in BHB fraction and temperature is 
relatively poor except at very old ages. 

We extend our analysis to show the same tracks (with 
the same color-coding) in two examples of color-spaces that 
combine UV and mid-IR data; FUV - NUV/NUV - [3.4] and 
FUV — NUV /NUV — r (Figure 11). These show a much wider 
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Figure 13. FSPS CSP templates depicting the evolution of FUV — NUV vs. NUV — [3.4] color-space with from ages 0.03-14. 1 Gyr. See Figure 10 for a full description. 


span in color on the x-axis by 5 mag, compared to the 0.7 mag 
range of the NUV— r/g — rand NUV— g/g — i plots (Figure 10). 

Additionally, there is a much cleaner separation in color 
for r, and somewhat cleaner separations for both BHB and 
temperature. The GALEX-WISE colors provide an expanded 
range in color-space for BHB fraction color-cuts, especially 
at high metallicity, and there is a marginal improvement for 
changes in A T . For example, comparing the top left plots in 
Figures 10 and 1 1, the different BHB fractions are separated on 
average by 1 mag in the GALEX-WISE, but only 0.2 mag in 
the GALEX- SDSS. A further advantage of NUV - [3.4] colors 
is the near all-sky coverage of both datasets. Because of these 
advantages, we adopted the FUV — NUV/NUV — [3.4] and 
FUV — NUV /NUV — r spaces as the primary color-space focus 
throughout the paper. 
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